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We consider the transverse Ising model in one dimension with nearest-neighbour interaction and 
calculate exactly the longitudinal spin-spin correlation for a class of excited states. These states 
are known to play an important role in the perturbative treatment of one-dimensional transverse 
Ising model with frustrated second- neighbour interaction. To calculate the correlation, we follow 
the earlier procedure of Wu, use Szego's theorem and also use Fisher-Hartwig conjecture. The 
result is that the correlation decays algebraically with distance (n) as l/\/n and is oscillatory or 
non-oscillatory depending on the magnitude of the transverse field. 
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I. INTRODUCTION 



Transverse Axial Next-Nearest Neighbour Ising 
(TANNNI) model has attracted much attention recently 
as a simple quantum Ising model with frustration. In 
one dimension, it is defined by the Hamiltonian [E 0, Q 
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N 



b j+2 



where sj and s| 



are the Pauli matrices a z and a x at 
the site j, J is the strength of the nearest neighbour 
ferromagnetic interaction, k is the frustration parameter 
and r is the transverse field. The phase diagram of the 
ground state of this model in the k-T plane shows in- 
teresting features and is a matter of much controversy 
and investigation. Recently, we have shown [4| that the 
ground state of the Hamiltonian Ha is related to a class 
of excited states of a much simpler Hamiltonian obtained 
by removing the frustration (i.e. setting k = 0), 
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This Hamiltonian defines what is called transverse Ising 
model (in one dimension) and is the first quantum Ising 
model to be solved exactly ([!, Q, see also 0, @])- 

All the eigenvalues and eigenfunctions of the Hamilto- 
nian H for transverse Ising model can be derived exactly 
P, H, IE 0] ■ The 2 N number of energy eigenvalues are 



E = 2T^ k A k 



(3) 



where A fc = ^(A 2 + 2Acosfc + 1), A = J/T, ^ may be 0, 
±1 and k runs over N/2 equispaced values in the interval 
to 7r. The ground state is obviously 
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corresponding to = — 1 for all k and the excited states 
correspond to the other distributions of values. The 
spin-spin correlation function in the longitudinal direc- 
tion 
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> - < s; 
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has also been calculated exactly 0, H[ for the Hamilto- 
nian H in the ground state; it decays exponentially with 
distance at A / 1 and algebraically (with exponent = 
1/4) at A = 1 (the critical point). As for the excited 
states, the correlation was first studied by McCoy 0]. 
The specific quantity calculated was the correlation at a 
temperature T > for the XY model (which in turn can 
be mapped onto the transverse Ising model). The cor- 
relation was found to decay exponentially with distance 
except when T = and A = 1 . 

The basic motivation of the paper is that (as mentioned 
above) the ground state of the Transverse Axial Next- 
Nearest Neighbour Ising Hamiltonian Ha in the region 
of interest (namely, k ~ 0.5, r ~ 0, sec rcf. [4]) is related 
to a particular set of excited states of the transverse Ising 
Hamiltonian H. This states have 



-1 for k — to <t> (unexcited) 
1 for k — (j> to 7r (excited) 



(5) 



It is thus important to study these states for transverse 
Ising model, especially, to calculate in these states, the 
longitudinal spin-spin correlation C z (n) for large n. In 
this paper we provide exact calculation for this quantity. 
Our result is that for these states, the correlation decays 
algebraically with exponent 1/2 for all values of A ^ 1, 
provided <f> is neither or 7r (some levels are excited, but 
not all). This communication does not cover the case of 
A = 1. 

Our results do not contradict the results obtained by 
McCoy, since at T > 0, not only the excited states de- 
scribed by Eq. (5) but all possible excited states con- 
tribute to the correlation. Thus, any state with energy 
E would provide a contribution by an amount propor- 
tional to exp(— E/kBT) (fee is the Boltzmann constant). 
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If the states with exponentially decaying correlation be 
dominant, the resultant correlation will not be algebraic. 

Instead of calculating C z {n) directly, we shall rather 
calculate the quantity 



a z (n) ~ 



(0) 



We shall see that for the excited states discussed in this 
paper (see Eq. (5)), 

lim a z (n) = 

n — >oo 

so that < s z > is also zero and C z (n) = a z (n). To eval- 
uate the quantity a z (n) for large n, we note that this 
quantity can be expressed as an n x n Toeplitz determi- 
nant dUi 
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where the elements are given by 
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In fact, Pfeuty [6j has given the expression for G m in the 
ground state and there this integral is from to tt. That 
expression has been modified here for the excited state 
following Lieb, Schultz and Mattis 9]. The standard 
prescription for studying the limiting (n — > oo) behaviour 
of a Toeplitz determinant is to apply Szego's Theorem in 
Fisher-Hartwig form [T7j|. [TTl [T^. Tl3l] . We shall also apply 
the same technique on the determinant given above. 

In Sees. II and III we treat the cases A < 1 and A > 1 
respectively and in Sec. IV we present some discussions. 



II. SPIN CORRELATIONS FOR A < 1 

We start this section with some analysis that is true 
both for A < 1 and A > 1 . The first step is to observe that 
the expression (8) for the elements G m of the determinant 
can be written as, 



1 f 2 * 

G m = — / dk exp(-imk) c(k) 
27r Jo 



(9) 



with the function c(k) (called the generating function) 
given by, 



w for < k < <fi 
c(k) = ^ —w for 4> < k < 2tt — <f> 
w for 2tt — 6 < k < 2tt 
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(ii) 



It should be mentioned that this generating function is 
similar to that for the classical two-dimensional Ising 
model |14| apart from the finite jumps at k = 4> an d 
2ir-<f). 

Till now, whenever we have come across a square root 
(Eqs. (3),(8),(11)) we have tacitly adopted the conven- 
tion (like usual computer language compilers) that we 
shall always take the numerically positive root. How- 
ever, since later (see Eqs. (16), (17) below) we shall need 
to integrate the logarithm of the generating function from 
to 2tt, let us choose, following the standard literature 
[TH , the branch of square root for which this function is 
positive in the limit k — > tt. Thus, we write the surd in 
Eq. (11) as 



A + cos k + i sin k 
v / (l + A 2 + 2Acos/c) 



(11a) 



and observe that we should take the negative root for 
A > 1 and positive root for A < 1. Let us call the gen- 
erating function evaluated with this new convention as 
c(k) which is related to the c(k) of Eq.(10) by 



-c(k) for A > 1 
c(k) for A < 1 



(12) 



This c(fc) when plugged in Eq. (9) leads to G m and 
the Toeplitz determinant det n (A) constructed from G m 
can be evaluated by applying Szego's theorem on c(fc). 
Obviously, this determinant is related to the det„(A) of 
Eq. (7) by 



det n (A) = 



(-l)"det„(A) for A > 1 
det„(A) for A < 1 
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The "index" of the generating function, defined as 



Ind c(Jfc) 
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[logc(27r) - logc(0)] 



is for A < 1 but +1 for A > 1. Since Szego's theorem 
can be directly applied only when the index is zero, our 
job is somewhat simpler for A < 1, compared to the same 
for A > 1. 

We shall be confined to the case of A < 1 henceforth. 
The generating function for this case is, 
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z(k) = < — w for <fi < k < 2n — <j) 
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Let us now investigate whether the conditions required 
(other than zero index) for Szego's theorem to be valid 
are satisfied here or not. For <fi — tt (ground state, all 
modes unexcited), one can make binomial expansion of 
the surd in w (Eq. (11)) and conclude that for any large 
positive integral value of m (with A < 1), 
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FIG. 1: 

Plot of the elements of the Toeplitz determinant for <j> — 0.8 
and A = 1.3. The oscillations themselves depend on <f> and A 
but the upper and the lower envelopes decay as 1/m for all 
A and 4>- The dotted lines are 0.635/m and 0.150/m. 



and 
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Since G m and G_ m both vanish exponentially with m, 
the conditions, 



E is, 



< oo 



and 
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are valid. Since c(k) #0 for all values of k, the Hirsh- 
mann conditions [16|, [T7| for the validity of Szego's the- 
orem is satisfied. (However, it is not clear whether the 
Devinatz condition [l6T| . namely, the continuity of 
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is true.) For </> ^ it, the asymptotic behaviour of | G m \ 
and | G-m | depends on A and <j> and is not easy to 
track analytically. Numerical study shows that both 
these quantities are oscillatory (Fig. [l}. The upper and 
lower envelopes of the oscillatory curve decay as 1/m for 
all A and <f> while the wavelength and amplitude of os- 
cillation depends on the specific values of A and <f>. It 
is thus difficult to analyse the validity of Hirshmann or 
Devinatz condition here. We do not proceed further with 
the question of validity of the conditions for applicability 
of Szego's theorem and rather choose to study the results 
obtained when we do apply the theorem anyway. 

To apply Szego's theorem, the next step is to calculate 
the quantities go and g± m defined as 



1 f 2 " 
•9o = 7T / dk mc ( fc ) 
2tt Jo 



(16) 



g±m = ^ r dke^ mk lnc(fc), 
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where m is any positive integer. These quantities can be 
easily evaluated from Eq. (13), but we must be careful 
to choose the principal branch for lnc(fc), satisfying 

—iir < lnc(fc) < iir. 

This criterion is satisfied if we express the logarithm as, 
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with the phase factor 



for < k < 4> 
-1 for 4> < k < 2tt - 4> 
-2 for 2vr - (j)< k < 2vr. 

By expanding the logarithm on the right hand side 
of Eq. (17a) in a series and performing a term-by-term 
integration, we now get 



and 
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.90 = 



(-A) m _^ cos(m(/>) 
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Next, we have to calculate the sum, 



a = E m -9™.9-m 
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(20) 
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Using the standard result 
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E— a m cos(ma) = -~ln(l- 2a cos a + a 2 ) (21) 
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(for < a < 2tt, a 2 < 1) one obtains, 



a=-ln(l-X 2 
4 v 



1 / 2sin0 

2 n Vl + 2Acos0 + A 2 
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The last sum is divergent. Following the conjecture of 
Fisher and Hartwig [ll|, [l^, [3i I3| > we replace the upper 
limit of the sum by C n to obtain 



E — = + 



(22) 



m — 1 



where is the Euler's constant (=0.5772). The constant 
C is independent of A, 4> and n. Its numerical value could 
not be found analytically. 

Szego's theorem now reads [T3, El, EH 

det„(A) = exp(ng + o), 
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which, on substitution of the expressions of go and <J leads 
to the final result 

C z n = det„ (A) = v =4= -= , 23 

VI + 2Acos</> + A 2 V™ 

where A = (l/^/C/2) exp(— 7^/2). We have numerically 
evaluated this constant to be 0.590 (which corresponds 
to C = 3.226) by computing determinants upto a size of 
200 x 200. The correlation function C z (n) decays alge- 
braically with distance without any oscillation. 



III. SPIN CORRELATIONS FOR A > 1 

We shall now evaluate the Toeplitz determinant of Eq. 
(7) for A > 1. Since the index of the generating function 
is not zero in this case (as mentioned above), Szego's 
theorem cannot be applied directly. One has to follow 
instead the method of Wu [20]. According to Eq. (8), 
the elements G m for some given value of A is related to 
the same for 1/A by, 



G TO (A) — G_ m +i(l/A) 



(24) 



Hence, the determinant for A = Ao > 1, can be written 
as. 
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Let us consider the Wiener-Hopf system of equations, 
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for k = 0, 1, 



■ n. By Cramer's rule, 
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This equation would give the value of the determinant 
for A > 1 in terms of the determinant for A < 1 provided 
we know x n for A < 1. 

To proceed further, we switch over to the convention 
that the root which ensures c(ir) > will be chosen and 
solve the Wiener-Hopf equations 



G/ s _ m (l/Ao) x m = Sk,o- 
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The solution for found by Wu [2lj is 
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with l(k) defined through |22j 
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Putting here the expressions of go and g± m (Eqs. (18), 
(19)) directly and using (21), we get 
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We rewrite this as 



= — dk cos(nfc) yl + A 2 + 2A cos k x 
1" Jo 
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— {cos(m0 + wfc) + cos(m0 — mfc)} 



exp 



and note that the formula (21) may be applied to 
^(1/m) cos(m(j)+mk) but not to XXl/ TO ) cos(m0 — mfc) 
since the former is convergent in the range of integration 
while the latter diverges at k = <j). Thus we have, 
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For large n, the term cos(nfc) oscillates rapidly and 
the non-trivial contribution to the integral comes only 
from the points where the rest of the integrand diverges. 
Hence, it suffices to integrate only over a narrow region 
(of width 2e <C 1/"-) around the point k — <p an d obtain 



1 + A 2 + 2Acos(</>) 



2tt 



sin(</>) 



cos(nx + «</ 



e 5 da; 



(34) 

where we have taken the smoothly varying part out of 
the integration and put x = k — <p. Following the Fisher- 
Hartwig conjecture, we now replace ( "heuristically" ) the 
upper limits in the sum S (see Eq. (33)) by C/n and 
note that for ~ 1, we must have xm <C 1 for all m 
in the range of integration. The sum S may then be 
approximated as 
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Substituting Eq. (22) in the first summation, we obtain, 
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As n — > oo, the portion within [•••] becomes the Dirac 
delta functon S(x) and the expression (34) for x n becomes 



e*> E v /l + A 2 + 2Acos(</>) 
0r sm(4>) 



cos(n4>) (37) 



The unknown constant £' (happily) does not appear in 
this final expression. Since x n = x n , Eqs. (12a), (23), 
and (27) now give our final expression for the correlation 
function for A > 1 : 
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(The constant A is defined below Eq. (23).) This cor- 
relation function decays algebraically with distance and 
also oscillates with distance with wave- vector (jr + <fi). 

Before we conclude this section, we point out that al- 
though Eq.(29) for derived by Wu [2l[ is applicable 
here, there is a major difference between our case and the 
one studied by Wu. In our case, x n does not vanish in the 
limit ti — ► oo and all we want to evaluate is the zero-th 
order term in n. But in Wu's treatment, Imin^^ x n is 
zero and the basic job was to evaluate the leading terms 
for large n. 



only for A = 1 that the correlation decays algebraically 
with distance. But for the type of excited states we 
study (see Eq. (5)), the algebraic behaviour is seen over 
a larger range, namely, for A < 1 as well as A > 1 and 
any <j> in the range < < 7r. 

(ii) For handling the divergent sums that arise in 
the calculation of correlation function, Fishcr-Hartwig 
prescription is very much useful for the excited states 
as it was for the ground state. Actually, the index for 
algebraic decay can be derived from the generating func- 
tion itself, without any further calculation. Fisher and 
Hartwig have shown that for generating function having 
multiple singularities, if arg[c(fc)] jumps by an amount 
2iriaj at the j-th singularity, then the determinant 
will decay algebraically with distance with an exponent 
J2j a2 j- In our case, the generating function of Eq. (10) 
has two jumps, at each of which exj is 1/2, and the 
index is 1/2, obeying this conclusion. Moreover, if we 
consider an excited state where (in contrast to Eq. (5)) 
the excited modes extend from k = <pi to fa, then (for 
4>2 7^ 7r) there will be four singularities in the generating 
function, at each of which aj is 1/2. The correlation will 
then decay as 1/n. 
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